Direct extraction of transversity and its accompanying T-odd distribution from 
the unpolarized and single-polarized Drell-Yan processes 
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The Drell-Yan (DY) processes with unpolarized colliding hadrons and with the single transver- 
sally polarized hadron are considered. The possibility of direct (without any model assumptions) 
extraction of both transversity and its accompanying T-odd parton distribution functions (PDF) 
is discussed. For DY processes measurements planned at GSI the preliminary estimations demon- 
strate that it is quite real to extract both transversity and its accompanying T-odd PDF in the PAX 
conditions. 
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The advantage of DY process for extraction of PDF, 
is that there is no need of any fragmentation func- 
tions. While the double transversely polarized DY 
process hIhI l + l~X allows to directly extract 
the transversity distributions (see ref.Q for review), 
in the single polarized DY H\H\ — > l + l~X the ac- 
cess to transversity is rather difficult since it enters 
the respective cross-section in the complex convolu- 
tion with another unknown T-odd PDF (see below). 
At the same time it is, certainly very desirable to man- 
age to get the transversity PDF from unpolarized and 
single-polarized DY processes as an alternative pos- 
sibility. Besides, T-odd PDF are very intriguing and 
interesting objects in themselves, so that it is very 
important to extract them too. 

The main goal of this paper is to investigate the pos- 
sibility to completely disentangle PDFs corresponding 
to the unpolarized and single-polarized DY processes. 

Let us first consider the results of ref. Q for both 
unpolarized and single-polarized DY processes. In 
that paper the Collins-Soper frame 1 is used ( see Fig. 
3 in ref. Q), where one deals with three angles 9, 4> 
and 4>Si- Two angles, and <f>, are common for both 
unpolarized and polarized DY processes. These are 
the polar and azimuthal angles of lepton pair. Third 
angle, (f>s 2 j does appear when hadron two is transver- 
sally polarized, and this is just the azimuthal angle of 
S 2 t measured with respect to lepton plane. 

We consider here the case of pure transverse po- 
larization of hadron two, so that we put Ai = and 
I Sir | = 1 (A2 = and |S 2 t| = 1 in our notation) in 
the respective equations of ref. Q (Eqs. (21) and (22) 
in ref. [2|) for unpolarized and single polarized cross- 
sections. Besides, taking into account only the dom- 
inating electromagnetic contributions and neglecting 
(just as in ref. (3) the "higher harmonic" term con- 
taining 30 dependence, one gets the following simpli- 
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1 See for detail of the respective kinematics. 



fied equations for the QPM unpolarized and single- 
polarized cross-sections : 
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Here h = qr/l^rl, h\ q {x,)/^) is the kx- 
dependent transversity distribution, while hj^Jx, k|,) 

and £j(x,ky) are fcr-dependent T-odd PDFs (see 
ref. 1] for review). The convolution product is de- 
fined 0] as 

T[f q f q ] = J d 2 k 1T d 2 k 2T S 2 (k 1T + k 2T - q T ) 

x [f q ( Xl , k 2 T )f q {x 2 , k 2 T ) + (1 «-> 2)] . (3) 

Let us first consider the purely unpolarized DY pro- 
cess. Notice that Eq. Q is very inconvenient in appli- 
cation because of the complicated qt and kr depen- 
dence entering Eq. via the convolution, Eq. ©. 



2 



To deal with Eq. Q the model 
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where Mn = 2.3 GeV , 
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1, ax = 1 GeF 2 and 



Mff is the hadron mass, was proposed in ref. Q. 
With a such assumption one then calculates 0, Q the 
coefficient k = j//2 at cos 20 dependent part of the 
ratio 
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which allows to explain 2 the anomalous cos 20 depen- 
dence 0,0 of the unpolarized DY cross-section. How- 
ever, the author of ref. |2( stresses that Eq. Q is just 
a "crude model". Besides, Eq. Q can not help us 
to extract the quantity hi from the unpolarized DY 
process. 

Thus, to avoid these problems, let us apply the qr 
weighting approach which was first proposed and ap- 
plied in refs. [6( and |jj with respect to a particular 
electron-positron annihilation process and in ref. [jj 
with respect to semi-inclusive DIS. To use the advan- 
tage of qr integration, one should extract from un- 
polarized DY process the properly integrated over 
ratio (c.f. Eq. ©) 
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that should be compared 3 with the equation (see refs 

HI) 

3 

R = ( 1 + A cos 6 + u sin 29 cos <j> 

167T 

+ (iV2)cos20sm 2 6>) [y = 2k, A ~ l,/x~0).(8) 

By virtue of Eq. Q), the coefficient k at cos 2(f) 
dependent part of R reads 



k = J d 2 q T [|q T | 2 /AfiM 2 ] 
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and, due to the properly chosen weight |qr| 2 , the inte- 
gration over qr leads 4 to the following simple equation 
for k: 



J2 q el(fi q (x 1 )h q (x 2 ) + (1^2)) ' ( } 



where the standard notation 



h^ n \x) = |d 2 k T ^) hi q (x,k$) (11) 

for the n-th moment of kr-dependent PDF is used. 
Thus, one can see that the numerator of k is factor- 
ized out in the simple product of the first moments 
of h^; distributions. This allows to directly extract 
these quantities from k which should be measured in 
unpolarized DY. This, in turn (see below), allows to 
directly extract the transversity distributions h\ from 
the single spin polarized DY. Notice that now there is 
no need in any model assumptions about kj- depen- 
dence of distributions. 

Let us now consider the single transversely polar- 
ized DY process HiH\ — > l + l~X and define the fol- 
lowing single-spin asymmetries (SSA) 

A Hf) = 

J dfld<j) S2 sin((f> ± 4>s 2 )[da(S2T) - da(S 2 T)} 

x^Jdnd<f>s 2 [d<T(S2T)+da(-S 2 T)]^ , (12) 

where the single-polarized cross-section is given by 
Eq. 10. It is clear that in the difference da(S2r) — 
da(— S 2 t) only the terms of Eq. $F2$ containing 
siri((f> — <ps 2 ) an d sin(</> + 0s 2 ) survive (and are multi- 
plied by two). Besides, the properly chosen 5 weights: 
sin(</> + cj>s 2 ) an d sin(0 — <ps 2 ), allow to separate the 
contributions containing h^ and f^j. PDF with the 
result 
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The asymmetries like Af given by Eqs. (|12|l . (|14|l 
and their application with respect to Sivers func- 
tion /^(a^kf,) = -(M/2|k T |)A^ HT (a:,k§,) extrac- 
tion from the data were considered in detail in refs. 



2 Notice that the large values of v cannot be explained by lead- 
ing and next-to-leading order perturbativc QCD corrections 
as well as by the high twists effects (see Q and references 
therein). 

3 Obtaining Eg. JSJ one sets || A = 1 and fj, = in the most 
general equation for R (Eq. (5) in ref. 0), which is justified 
Q by the expectation from next-to-leading order QCD and 
the data (refs 00) in the Collins-Soper frame. 



4 The normalization condition J d 2 ]ij'fi q (x, ky,) = fi q (x) is 
used (see, for example ref. Q). 

5 The analogous weighting procedure was applied ;9!| in the case 
of transversely polarized SIDIS by the HERMES collabora- 
tion. 
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[Tol[T]| . so that we concentrate here on the asymmetry 
A h given by Eqs. (JT^J) and lfl3l . 

Notice that asymmetry Ah given by Eqs. (fl"2*)l . (flTTft 

is inconvenient in application because of the compli- 
cated qr and kx dependence entering the convolution. 



So, we again apply the qx integration method 
(see also its application for the SIDIS processes in ref. 
9] and for the Sivcrs PDF extraction from the single 
polarized DY in ref. 0] ): 
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so that one easily gets 
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Thus, one can see that Ah is also factorized in the 
simple product of hf^ and hi. 

Among variety of DY processes, DY processes with 



antiproton (pp — ► l + l~X, pp^ — > l + l~X, p^p^ — > 
l + l~X) have essential advantage because the charge 
conjugation symmetry can be applied. Indeed, due 
to charge conjugation, antiquark PDF from the an- 
tiproton are equal to the respective quark PDF from 
the proton. Thus, Eqs. I|10[l. i|16[l in the case of pp 
collisions are rewritten as 
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where now all PDF refer to protons. Neglecting 
squared antiquark and strange quark PDF contribu- 
tions to proton and taking into account the quark 
charges and u quark dominance at large 6 x, Eqs. i|17f) 
and i|18fl are essentially given by 



k(xi,x 2 ) 
and 
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One can see that the system of Eqs. (|T§)) and (|2*U|l has 
very simple and convenient form in application. Mea- 
suring the quantity k in unpolarized DY (Eqs. ©, 



0) and using Eq. I|19|) one can obtain the quantity 
h^ 1 '. Then, measuring SSA, Eq. (fTST) . and using 

the obtained quantity hf^ , one can eventually ex- 
tract the transversity distribution h iu using Eq. lf2*U|) . 
Let us stress once again that now there is no need in 
any model assumptions about kx dependence of 
distributions. 



In order to obtain squares of and fi u in Eqs. 

(|19l) and (|2U[) . one should consider them at the points 7 
x\ = x 2 = x (i.e., xf = xi — x 2 — 0), so that 



iX(i)/ s , , , . k(x,x) 
hiu \ x ) = hu{x)\ — - — : 
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The large x values is the peculiarity of the pp experiments 
planned at GSI — see ref. hj 



7 The different points xp = 
value at fixed s = iiio Q 2 



can be reached changing Q 
= ^r>2 



4 



and 

h lu (x) = -4V2 A *h*L f lu (x), (22) 

yk(x,x) 

To estimate the possibility of and h\ u mea- 

surement, the special simulation of DY events with 
the PAX kinematics [12j are performed. The proton- 
antiproton collisions are simulated with PYTHIA 
event generator Two samples are prepared: for 
the collider mode (15 GeV antiproton beam colliding 
on the 3.5 GeV proton beam) and for fixed target 
mode (22 GeV antiproton beam colliding on an in- 
ternal hydrogen target). Each sample contains about 
100 K pure Drell-Yan events. Notice, that this is just 
the statistics planned to be achieved by PAX. Indeed 
(see ref. H2), ^he sample for collider mode corre- 
sponds to about one year of data-taking with a cross- 
section of 40 mb and a luminosity of 2 x 10 30 cm _2 s _1 . 
For fixed target mode it can takes about three months 
with a cross-section of 30 mb and a luminosity of 
about 10 31 cm _2 s _1 . 

Unfortunately, the original PYTHIA generator we 
deal with does not reproduce the corresponding to DY 
experiments 0, nontrivial qx and x dependencies of 
the quantity v entering Eq. (J8J. So, to estimate the 

possibility of and h\ u measurement, one should 

properly introduce these dependencies in accordance 
with the existing experimental data. To this end we 
apply the commonly used Monte-Carlo method based 
on weighting of the kinematical events. To apply the 
weighting procedure in our case, we just ascribe to 
each event the weight w — R which, in accordance 
with the data is given by Eq. JSJ), where A ~ 1, 

fi ~ and v has nontrivial and x dependencies. 
The qr dependence of v is taken from refs. |2j, |3j - 
Eq. (49) in ref. @ and Eq. (21) in ref. 0, and this 
qt dependence properly fits the existing experimen- 
tal data 0, 0| . However, in refs. 0, y| (where the 
simplified Boer's model is applied) there is no (im- 
portant and corresponding to DY experiments 0,0) 
cc-dependence of v at all, so that we take this depen- 
dence from ref. Q. 

To check the validity of the angular distribution 
analysis of the weighted events we reconstruct the qr 
and x\ dependencies of v. The results are shown in 
Figs. 1, 2. One can see a good agreement 8 between 
input (solid lines) and reconstructed (points with er- 
ror bars) values. 

Thus, applying the above described weighting pro- 
cedure, our simulations reproduce the nontrivial an- 
gular dependence of R with qx- and x-dependent v. 
These dependencies are in accordance with the respec- 
tive input dependencies obtained in experiments on 



As an additional check of our analysis validity, we reproduce 
the input zero value of u. 
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FIG. 1: Reconstructed from simulations (fixed target mode) 
quantities fi and v versus qx in comparison with the input 
(corresponding to experimental data) dependencies (solid 
lines). 
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FIG. 2: Reconstructed from simulations (fixed target mode) 
quantities /i and v versus x\ in comparison with the input 
(corresponding to experimental data) dependencies (solid 
lines). 

DY 0,0- Now it is straightforward to reconstruct the 
qT-weighted quantity R (Eq. EJ) and, consequently, k 
(Eq. [7J|. The results are shown in Fig. |3| The values 
of k at averaged Q 2 for both modes are found to be 
1.2±0.2 for collider mode and 1.0±0.2 for fixed target 
mode. 

The quantity is reconstructed from the ob- 

tained values of k using Eq. lj?T|) with xp = ± 0.04. 
The results are shown in Fig. 0] The obtained 
magnitudes of hf£ are in accordance (in order of 
value) with the respective magnitudes obtained with 
the model (0} for h^ u (x, k^). Indeed, for example 
for the collider mode {Qa Verage — 9 GeV 2 , so that 
x\ ~ X2 — 0.2 at the point xp ~ 0) the results from 
the simulations and from the model (@J are h^ 1 ' ~ 1 
and fcy^ — 0.5, respectively. 

Using the obtained magnitudes of we estimate 
the expected SSA given by Eq. (J20J. The results are 
shown in Figs. 5 and 6. For estimation of h\ u entering 
SSA together with (see Eq. (O) we follow the 
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FIG. 3: k versus x\ at xf — 0. Data is obtained with 
MC simulations for collider (closed circles) and for fixed 
target mode (open circles). For better visibility (to avoid 
overlapping) the points for collider (fixed target) mode are 
shifted 0.01 to the left (right) along the x-axis. 
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FIG. 4: h llt versus x\ at xf — 0. Data is obtained with 
MC simulations for collider (closed circles) and for fixed 
target mode (open circles). For better visibility (to avoid 
overlapping) the points for collider (fixed target) mode are 
shifted 0.01 to the left (right) along the x-axis. 



procedure of ref. |l4j and use (rather crude) "evolu- 
tion model" > where there is no any estimations 
of uncertainties. That is why in (purely qualitative) 
figures 5 and 6 we present the solid curves instead 
of points with error bars. To obtain these curves we 
reproduce x-dependence of hf^ in the considered re- 
gion, using the Boer's model, Eq. (@J, properly nu- 
merically corrected in accordance with the simulation 
results. 

To estimate the measurability of the quantities we 
deal with, it is relevant to estimate the upper bounds 
on hi, hj) and then on k and Ah- Obtaining 
hf^ and hi u one deals with Eqs. (jTT?)l . lj2U|) ap- 
plied at the points xi ~ X2 — \J Q 2 / 's, so that we 
perform the estimation of the upper bounds on k 
and A/j at the points Xp ~ corresponding to the 
average Q 2 values for both collider and fixed target 



modes. The maximally allowed value of h 1 can 
be found operating just as it was done with respect 

to the quantity (first moment of the Sivers 

function) in refillOj. To this end we first apply 
the inequality 9 [l5j(|kT | /M)hj; (x,\s%)< fi(x,\s%). 
Then, using the estimation (see ref. [T(j and refer- 
ences therein) (kr) — 0.8 GeV one easily gets the up- 
per bound on hf^: /i^ < 0.4/i„(x). On the other 
hand, maximally allowed value of hi u can be found 
using the Soffer inequality |fti u | < (fi u + gi u )/2. 
For the PAX kinematics s = AZGeV 2 , Q 2 average ^ 
5 GeV 2 for the fixed target mode and s = 215 GeV 2 , 
Qaverage — 9 GeV 2 for the collider mode. Thus, at 
the point xp — we deal with, xi ~ X2 — 0.3 and 
x\ ~ X2 — 0.2 for the fixed target and collider modes, 
respectively. Then, the inequalities on hi u and 
give 10 h lu(max) ~ 1.5 {fin = 1.9, g lu = 1.0 ) and 

h lu(max) — °' 8 for fixed tar S et mode wnile hiu (max) — 

2.3 (/ lu = 3.1, g lu = 1.5) and hf^ max) ~ 1.2 for 
collider mode. Using these estimations of hi u ( max ) 
and h^^ max ^ in Eqs. ljl"§|) . (f^Oft it is straightforward 

to obtain the maximally allowed values of k and Ah'. 
k(max) ^1.4 and \A h{max )\ ~ 0.17 for fixed target 
mode while k (max) ~ 1.2 and \A h(max) \ ~ 0.14 for 
collider mode. One can see that obtained estimations 
of upper bounds on hf£ , k and Ah are in accordance 
with the results presented by Figs. 3-6. 

Looking at the (preliminary) estimations presented 
by Figs. 3 and 4, one can conclude that the quantities 
k and are presumably measurable in most of 

the considered rr-region. At the same time, looking 
at Figs. 5 and 6 one can see that for both modes 
SSA Ah is estimated to be about 6-8%. On the other 
hand, as it was argued in Ref. 0] (see section "Single 
Spin asymmetries and Sivers Function", p. 25), the 

studied in ref. [13 SSA A^ of order 5-10% 

can be measured by PAX. It is obvious that studied in 
this paper SSA Ah, weighted with sin(0 + 0s) and the 
same weight qt/Mjv, is absolutely analogous to SSA 

sin(0— (fis) 2if sin((p—<ps ) 

A UT N , so that it is clear that if A UT N 

of 5-10% is measurable, then Ah of 6-8% is measurable 
too. 

Thus, it is shown that it is possible to directly ex- 
tract the transversity and its accompanying T-odd 
PDF from the unpolarized and single polarized DY 
processes with antiproton participation. It is of im- 
portance that there is no need in any model assump- 



9 This inequality i s di rectly obtained by relaxing the bound 
Eq. (16) in ref. 1158 (eliminating the unknown distribution 
in that bound). 

Performing the se estimations we use GRSV2000LO 
parametrization Il7l for g\ u and GRV98LO parametrization 
El for f Xu . 
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FIG. 5: SSA given by Eq. 12(A) versus xf for collider mode 
for three values of Q 2 : 50 GeV 2 (lower curve), 25 GeV 2 
(middle curve) and 9 GeV 2 (upper curve). 



tions about kx dependence of hj~ . One can directly ex- 
tract both hi and first moment of /ip from the single- 
polarized and unpolarized DY processes, since these 
quantities enter the measured k and SSA Ah in the 
form of simple product instead of complex convolu- 
tion. The preliminary estimations for PAX kinemat- 
ics show the possibility to measure both k and SSA 
Ah and then to extract the quantities and hi. 

Certainly, the estimations of k and Ah magnitudes 
obtained it this paper are very preliminary and show 
just the order of values of these quantities. For more 
precise estimations one needs the Monte-Carlo gener- 
ator more suitable for DY processes studies (see, for 
example, ref. Q ) than PYTHIA generator which we 
used (with the proper weighting of events) here. 
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FIG. 6: SSA given by Eq. \2(A) versus xf for fixed tar- 
get mode for three values of Q 2 : 16 GeV 2 (lower curve), 
9 GeV 2 (middle curve) and 4 GeV 2 (upper curve). 



Notice, that it is straightforward to properly modify 
the procedure discussed in this paper to DY processes: 
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be study [Hj in the COMPASS experiment at CERN. 
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